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ABSTRACT 

This paper^describes a student project in digital 
simulation techniq'ues that is part of a graduate systems analysis 
course entitled Biosimulation. The students chose different 
simulation techniques to 'solve a. problem related to the neuron model. 
(HLH) 
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Introduction ' ^ 

Among the graduate offerings of the Institute for Sensory Research in the 
area of systems analysis is a course in Biosimulation. An integral part of 
the course is a student project devoted to the application of digital simula- 
tion techniques. Difficulties arise in the choice of these projects because 
of the wide diversity of student backgrounds and because of different levels 
of computer skill. I would like to share with you my experience in attempting 
to-x>vercome some of these difficulties. 

As a student project we chose a relatively simple^ but realistic, problem 
having no known analytic solution. Concentration. would be on the simulation 
rather than on difficulties inherent in a complex problem. All students 
solved the same problem but were allowed to choose a digital simulation 
technique based on their individual backgrounds, level of computer expertise, 
availability of computer resources, and their personal interest in the prob- 
lem. Class format was essentially a seminar with e^ch student giving an 
oral presentation midway in the project and at completion. The mid-project 
presentations provided an opportunity for students to compare solution 
techniques and implementations, and to discuss real or imagined obstacles. 
Final presentations were aimed at comparing computed Results as well as 
discussing the effectiveness of individual simulations in terms of direct 
cost, accuracy, and programming difficulty. 

The Problem 

The problem chosen involves firing dynamics of the simple neuron model 
shown in Fig. 1. Specifically, the simulation was to determine characteris- 
tics of the spike train produced when summated excitatory post synaptic 
potentials (EPSPs) interact with a dynamic threshold. EPSPs, each with an 

.3 



2 



exponential tail of known time constant, occur randomly with inter-event 
times controlled by a Poisson process. An output spike is generated when 
the summed EPSPs exceed threshold. After each spike an absolute refractory 
period (ARP) is entered during which time the threshold is infinite. 
Following the ARP, threshold recovers as K/t toward its resting value. 

Implementations 

Figure 2 summarizes the diversity of stud^t backgrounds as well as 
the various co^pputers and program type3 chosen. The large computers, 
IBM 370 and PDP 10, are part of the Syracuse University Computing Center and 
were used in both the batch and terminal modes. The LINC 8 computer is 
located within the Institute for Sensory Research, where many of the students 
are graduate assistants, and was available at no cost. 

Algorithms of three main types were chosen. They included a straight- 
. forward calculation that proceeded point-by-point in time; an interesting 
time calculation that skipped over portions of time where threshold crossings 
could not occur; and a rather novel table look-up program based on the shape 
invar iance of a sum of exponentials of equal time-constants. More adventure- 
some students combined various features of these algorithms. Simplified 
flow charts for the "Modified Interesting Time" and/^able Lookup" programs 
are shown in Figs. 3 and 4 respectively. 

/ In the Modified Interesting Time program, a calculation ig made at 
adjacent sample points to determine geometrical conditions. Five conditions 
are recognized: onset crossing, diverging, converging, intersample crossing, 
•and possible tangency. In the latter two situations, a smaller sampling 
interval is chosen to more closely determine the time of threshold crossing. 
Thus the program dynamically selects an appropriate sampling interval. When 
* 

4 



ERIC 



conditions exclude the possibility ot a threshold crossing, the program 

advances time to the next interesting event, namely the onset of a new EPSP. 

The Table Lookup program of Fig. 4 takes advantage of the fact that 

the sum of exponentials , each with the same time-constant, is an exponential 

having the same time-constant but a different scale factor. The 'program 

-mm 

consists of two permanent tables. One stores valties of the exponential 
function at equally spaced points in time; the second stores values of the 
absolute refractory period and the K/t recovery function using the same time 
increment. Pointers are set to the head of each table and the tables com- 
pared entry by entry. If a crossing occurs (THD entry » EPSP entry times 
EPSP multiplier), the threshold pointer is returned to the head of its table, 
an output event is generated, and the coiiq)arison between tables is continued. 
When a new EPSP arrives, the EPSP multiplier is updated and, the EPSP pointer 
is returned to the head of its table. Coiiq)arisbn between tables continues as 
before. This Table Lookup program is essentially a point-by-point program 
but with a particularly fast and simple computational algorithm. 

Simulation Outputs and Results 

Output formats varied among students, and ran the ^gamut from exhaustive 
tabular values of threshold and generator potential at each sample point to 
summaries based on spike generation times. Ail example of the latter, type is 
shown in Fig. 5. This output is from a run made with the Modified Interesting 
Time program. The number of EPSPs to be processed and the average rate of 
occurrence are entered as parameters. Each line of the printout corresponds 
to the generation oi a spike. Each X represents an EPSP that has occurred 
since the last spike. The number at the end of the line is the inter-spike 
interval, the output report concludes by giving the total number of EPSPs, 
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the total number of spikes generated, the time period simulated, the average 
EPSP rate, and the average firing rate. \^ 

A completely different output report was generated using the CSMP 

(Continuous System Modelling Program) program. This Is a llne-prlnter 

u \ 
gifaphlcal display of hoth the generator potential and threshold as a function 

of time. Figure 6 shows a short segment of this report covering about 10 
milliseconds of real time. The left column is the generator potential as a 
function of t^me; the right column is the threshold. The high flat portions 
occur during the ARP. Paging format in the line pr4.nter creates the blank 
portion across the center. 

A very exciting and satisfying aspect/ o"?s^he project occurred when in- 
dividual results were pooled and plotted./ Figure^? is a pooled plot of 
average firing rate as a function of average EPSP rate. Results for all 
seven students are plotted. The only major deviation is for the .data 
represented by open diamonds. The fall-^ff at high EPSP rates is due to 
sampling errors. 

Conclusions , 

This student project restulted in a number of accoiE?)llshments, First, 
it acted as a vehicle to uncover many aspects of digital simulation. It 
was an introduction to practical computing for some; an extension of computing 
knowledge for others. It allowed the comparison of a number of high-level 
computing languages and programming techniques. It provided a practical 
seminar in oral presentation and in program documentation. Finally, the 
actual results obtained are a step in the solution of an Important problem 
in neural dynamics related to on-going research at the Institute. 
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The project generated an unusual amount of student enthusiasm and 
provided experience in. a significant and non-trivial sinrulation study. Each 
student, regardless of background, was able to achieve easily comparable 
minimal results, yet all were challenged at their own level of competence. 
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